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Abstract 

It is argued that the very first signatures of the approach to the 
black disk asymptotical limit in hadron collisions may be observed in 
the differential cross section of elastic scattering. The exponentially 
decreasing with the angle (or regime beyond the diffraction peak 
will become replaced by an oscillatory behavior or by the power-like 
falloff. Some estimates of energies where this can happen are pre- 
sented. 

PACS: 13.85.Dz 

Hadrons behave as extended objects in high-energy collisions. The ex- 
treme model widely discussed for many years is the asymptotical picture of 
two absolutely black Lorentz-contracted disks with complete absorption at 
any impact parameter. In the impact parameter formalism the total cross 
section is given by the expression 

<r t (s) = 2 J d 2 b(l-e- n{s ' h) ). (1) 

The integration is over the two-dimensional space of the impact parameter 
b. Q(s, b) denotes the opaqueness (or blackness) of disks colliding with 
the c.m.s. energy yfs and the impact parameter b. For absolutely black 
(fi(s, b) — > oo, b < R; Q(s, b) = 0, b > R) and logarithmically expanding 
(R = Ro In s, i? =cons t ) disks with radii R one gets from ([I]) 

a t (s) = 2nR 2 + 0(\ns). (2) 

The elastic and inelastic processes should contribute on equal footing 

-^ = ^\ = \tO{\^s). (3) 
a t {s) a t (s) 2 

The width of the diffraction peak I? -1 (s) should shrink because its slope 
increases as 

B(s) = ^ + 0(\ns). (4) 



The forward ratio of the real part to the imaginary part of the elastic 
scattering amplitude p must vanish asymptotically as 

p (s) = £- s +O(\n~ 2 s). (5) 

None of these asymptotical predictions were yet observed in experiment. 

Total cross sections increase with energy albeit all fits taking into account 
the latest LHC data pQ prefer the power-like behavior. The preasymptotic 
origin may be ascribed to this effect, especially in view of some theoretical 
guesses [2] proposed long time ago. 

Inelastic processes are about 3 to 4 times more abundant than elastic 
ones. Even though the elastic cross section increased from 18% of the total 
cross section at ISR energies to 25% at LHC energies, it is much less than 
50% shown above. 

The diffraction peak shrinks about twice from energies about y/s « 6 
GeV where B ns 10 GeV~ 2 to the LHC energy where B « 20 GeV~ 2 . At 
ISR energies the slope B(s) increased logarithmically. To account for LHC 
data one needs a stronger than simple logarithmical dependence. The terms 
proportional to In 2 s are usually added in phenomenological fits. Even then 
the predictions [3j H] are not completely satisfactory. In connection with 
the power-like preasymptotical behavior of at one could expect at present 
energies the faster than logarithmical shrinkage of the diffraction peak as 
well. 

As a function of energy, the ratio po increases from below, crosses zero and 
becomes positive at high energies. This is a general tendency for collisions 
of any initial particles. For pp-scattering, the prediction of (jSJ) with values 
of s scaled by 1 GeV is still somewhat higher (about 0.177) than estimates 
from dispersion relations (0.14 in Refs [5], E]) even at 7 TeV while strongly 
overshoots them at ISR where tt/ Ins ~ 0.37. No logarithmic decrease is seen 
in these predictions which, however, depend on behavior of the total cross 
section at higher energies. Moreover, the value 0.14 can only be reached 
according to (jSJ) at the energy 75 TeV. Probably, at energies higher than 75 
TeV the first signs of approach to asymptotics will become visible. No data 
about po at the LHC energies exist yet. 

Thus we can only hope that we are now in the preasymptotical regime, 
at best. 

Surely, there is another more realistic possibility that the black disk model 
is too extreme and the gray fringe always exists. It opens the way to the nu- 



merous speculations with many new parameters about the particle shape and 
opacity (see, e.g., [H El El HOI HH H21 HH H31 HS1 HS1 HZ] which, nevertheless, 
are not as extreme as the black disk model. 

However, let us study just the black disk model. There exists another 
parameter which can shed some light on the situation. This is the ratio 



(6) 

which should be equal to 1/2 for black disks. In the Table 1 we show how 
this ratio evolves with energy. All entries there except the last two are taken 
from Ref. [7] with simple recalculation Z = 1/AY. The data at Tevatron and 
LHC energies are taken from Refs [T8~| [1] . All results are for pp-scattering 
except those at 546 and 1800 GeV which are for pp processes. The accuracy 
of the numbers listed in the Table 1 can be very approximately estimated as 
±0.1 from known error bars for the total cross sections and the slopes. 

Table 1 

v/£, GeV— 2.70 4.11 4.74 6.27 6.98 7.62 13.8 62.5 546 1800 7000 
Z - 0.64 1.02 1.09 1.26 1.31 1.34 1.45 1.50 1.20 1.08 1.00 

We see that Z increases from low energy region to ISR energies and then 
drops down approaching 1 at LHC. Our interest in this ratio is related with 
the fact that the differential cross section of elastic scattering beyond the 
diffraction peak in the so-called Orear region is very sensitive to its value. 

This region can be well fitted at different energies from rather low ones [19] 
to TOTEM data at LHC [20] by the very general expressions describing the 
cross section falloff with the dip-bump substructure derived in Refs [211 122] . 
It was shown there that the unitarity relation outside the diffraction peak can 
be reduced to the integral equation for the imaginary part of the amplitude. 
No assumptions have been made other than the validity of experimental data 
in the diffraction peak. Using the shape of the amplitude in the diffraction 
cone 

A d (p,9)=Atp 2 a t e- B ^ 2 / 2 (l~ iPd ), (7) 
one gets the unitarity equation written as 

lmA(p, 0) = V °* dQ x e~ B ^ e - e ^l\\ + PdPl )lmA(p, 6,) + F(p, 6). 

A<k\/2ixB J-oo 

(8) 



The amplitude is denoted by A. p« a/s/2 and 9 denote the momentum and 
the scattering angle in the center of mass system, p^'s take into account the 
real parts at the corresponding angles in the diffraction cone (i = d) and at 
larger angles (i = I). The integral term represents the two-particle interme- 
diate states of the incoming particles. The function F(p,8), called following 
Ref. [23] as the overlap function, represents the shadowing contribution of 
the inelastic processes to the elastic scattering amplitude. 

Its analytic solution was obtained in a model-independent way with the 
assumptions that the role of the overlap function F(p, 9) is negligible outside 
the diffraction cond3 and the real parts may be replaced by their average 
values in the diffraction peak pd and outside it pi, correspondingly. Let us 
stress once more that the Gaussian shape of the amplitude has been only used 
at rather small angles in accordance with experimental data. If it is ascribed 
to the Pomeron exchange, one is tempted to consider this solution as an effect 
of shadowing corrections due to multi- Pomeron graphs. Such corrections near 
the diffraction cone were considered, e.g., in Refs. [25, 26J. Even though 
their conclusions differ, the qualitative tendency of peak suppression is valid. 
The definite model summing multi- Pomeron graphs was considered in [27] . 
However, the oscillations over the exponentially damped regime happened 
to be too strong. The non-linear unitarity equation was solved by iterations 
assuming the Gaussian shape of the overlap function [22] but did not lead 
to quantitative agreement with experiment. That is ascribed to the ad hoc 
chosen shape of the overlap function which was recognized as unsatisfactory 
in earlier works (e.g., see [281 12H1 [301 EI])- The solution of Eq. (jSJ) does not 
depend on model assumptions. 

Here we show only the leading term of the solution describing the expo- 
nential falloff with the angle 9 (or ^/\t\) within the Orear regime which is 
most important for our purpose: 



lmA(p,9) = C (p)exp -J2B\n— p9 . (9) 

\ V ^ + pdpi j 

The solution predicts the dependence on p9 m y\t\ but not the dependence 

1 The results of the papers [2Ql [24] give strong support to this assumption. The overlap 
function has been calculated there directly from experimental data subtracting the elastic 
contribution I2 from the lefthand side of the unitarity equation wthout model assumptions. 
It is extremely small outside the diffraction peak, and becomes even smaller at the LHC 
enargy compared to the lower ones. 



on the collision energy. Its validity may be checked by the direct substitution 
in Eq. (jSJ). The values of the parameter Z discussed above are crucial 
for the slope of the differential cross section in the Orear region. Another 
important parameter is the product of the average values of ratios of real 
to imaginary parts of the amplitude in the diffraction cone and in the Orear 
region. While the former ratio can be estimated from the experimental results 
on interference of Coulomb and nuclear amplitudes supported by calculations 
with the help of the dispersion relations, the latter one was never estimated 
from experimental data. Any information about it becomes valuable. That 
is why it is important that from the fit [20] at 7 TeV it became possible for 
the first time to get its estimation from the experimentally measured slope 
in the Orear region. It happened to be surprisingly large in the absolute 
value (pi ~ —2). The high sensitivity of the expression (Q to this parameter 
is determined by the value of Z being practically equal to 1 as seen in the 
Table 1 and, correspondingly, of InZ very near 0. 

Once the expression for the imaginary part of the amplitude (Q outside 
the diffraction peak is known, it is possible to calculate the real part in the 
same region assuming the geometrical scaling [32] is valid. This has been 
done in [33J. It is shown that the real part changes its sign in accordance 
with the general theorem proven in [31] and becomes negative in the Orear 
region. Its absolute values happen, however, to be much smaller than those 
required by the fit of TOTEM data at 7 TeV. This problem has not been 
resolved yet. 

Now, the question arises about what happens if the black disk limit of 
Z — 1/2 will be approached at higher energies. If the exponential slope 
persists, the only solution would be to use even higher absolute values of pi 
(but negative!). However, this would show that there is no approach to black 
disks at all. This limit asks for p d oc In -1 s and, consequently, to pi oc In s 
that seems too exotic possibility. If the product pdpi tends to zero that 
looks more realistic, then the behavior of the differential cross section in the 
Orear region should become cardinally changed. In place of the exponentially 
decreasing falloff, the oscillatory behavior should appear. Now, the solution 
of the equation (|S]) looks as 



lmA(p, 6) = C (p) cos 
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(10) 



The imaginary part of the amplitude dominates and determines the shape of 



the differential cross section 
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(11) 



The oscillatory behavior with periodically placed zeros is clearly seen. 
In the black disk limit Z = 1/2; p t = 0; B = R 2 /A. Thus 

lmA(p, 6) = C (p) cos^O^ 2 ln2|t|) = C (p) cos(^a t ln2 /4:iip9). (12) 

The distance between oscillation zeros becomes very small and asymptotically 
tends to zero 



It is hard to predict at what energies this transition from one regime to 
another will happen. The only guess I can propose is to rely on possible 
approach of the values of pd ~ p Q evaluated from dispersion relations to the 
black disk limit (JSJ). As discussed above, it can happen at energies about 75 



In some way, the conclusion about the undamped oscillations looks con- 
tradictory. Hardly the differential cross section will not decrease with \t\ at 
all. It is more probable that the factor Cq/s 2 is so strong that the regime 
( ITTj) dies out very fast with energy increase. Then it becomes unobservable 
and the place of the Orear region is filled by the hard parton scattering. The 
exponential with |t| 1//2 falloff will be replaced by the power-like decrease. 

Anyway we conclude that the approach to the black disk limit will induce 
drastic changes in the behavior of the differential cross section in the Orear 
region replacing the exponential falloff either by undamped oscillations in 
some energy interval or, that is more probable, by power-like behavior due 
to the hard parton scattering. It is possible that this feature will become 
the very first one noticeable if collision partners remind the black disks at 
extremely high energies. 
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